High count rate y-ray spectroscopy with LaBr3:Ce scintillation detectors 



B. L6her* Ab , D. Savran a - b - c , E. Fiori ab , M. Miklavec c , N. Pietralla d , M. Vencelf 

"ExtreMe Matter Institute EMM1 and Research Division, GSI Helmholtzzentrum fur Schwerionenforschung, Planckstr. I, 64291 Darmstadt, 

Germany 

b Frankfurt Institute for Advanced Studies FIAS, Ruth-Moufang-Str. I, 60438 Frankfurt am Main, Germany 
c Innovation Centre for Advanced Sensors and Sensor Systems, INCAS 3 , Assen, The Netherlands 
" Jozef Stefan Institute, Jamova cesta 39, 1000 Ljubljana, Slovenia 
' Institut fiir Kernphysik, TU Darmstadt, Schlossgartenstr. 9, 64289 Darmstadt, Germany 



Abstract 

The applicability of LaBr3:Ce detectors for high count rate y-ray spectroscopy is investigated. A 3"x3" LaBr3:Ce 
detector is used in a test setup with radioactive sources to study the dependence of energy resolution and photo peak 
efficiency on the overall count rate in the detector. Digitized traces were recorded using a 500 MHz FADC and 
analysed with digital signal processing methods. In addition to standard techniques a pile-up correction method is 
applied to the data in order to further improve the high-rate capabilities and to reduce the losses in efficiency due to 
signal pile-up. It is shown, that y-ray spectroscopy can be performed with high resolution at count rates even above 
1 MHz and that the performance can be enhanced in the region between 500 kHz and 10 MHz by using pile-up 
correction techniques. 
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1. Introduction 

In modern nuclear physics experiments y-ray spec- 
troscopy is one of the major tools to investigate the 
structure of stable and unstable nuclei. The detec- 
tion and spectroscopy of the emitted photons in the de- 
excitation of particle-bound excited states is a manda- 
tory ingredient in many studies of nuclear reactions. 
Consequently y-ray detectors play an important role in 
the configuration of experimental setups. 

Different types of detectors are available. The op- 
timum choice depends on the requirements of the y- 
ray spectroscopy experiments. High-Purity Germanium 
(HPGe) detectors are often used for highest resolution, 
while several scintillator based detectors provide better 
timing properties (at the cost of reduced energy resolu- 
tion) or are more suited for high acceptance spectrome- 
ters due to their much lower costs. 

In many experiments, the emitted y rays in the reac- 
tion of interest are accompanied by a rather large y-ray 
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background close to the target position. This radiation 
originates from photons produced in background reac- 
tions such as bremsstrahlung by fast secondary elec- 
trons or reactions within the beam dump or other ma- 
terial along the beamline. Since it is difficult (or even 
unwanted) to suppress this intense background one of- 
ten has to deal with high count rates. To avoid signal 
distortion, the maximum count rate capability of each 
detector must not be exceeded. As a consequence the 
rate in the photon detectors often limits the beam inten- 
sity and therefore the reaction rate. 

The recently developed scintillator mate- 
rial lanthanum bromide (LaBr3:Ce, tradename 

TM r~L 

BrilLanCe 380 by Saint Gobain [1]) is well suited 
to push the limitations posed by this problem. Its 
excellent timing properties (down to lOOps in optimal 
conditions 12j-|4[]), temperature stability 15|,|6[] and high 
energy resolution for a scintillator (3% at 662 keV 
|7,|8J]), which is dominated by statistical contributions 
091], have made LaBr3:Ce the material of choice for 
many nuclear ph ysic s experiments includin g y -ray 
spectroscopy |l0l yjj], medical imaging B12l 1 1 311 or 
industrial applications lfl4l [l5ll . The high light yield 
(165% of that of Nal) and the very short scintillation 
light decay constant (between 20 and 30 ns J31 Si) 
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Figure 1 : A schematic illustration of the experimental setup. The 
weaker I37 Cs source was kept at a fixed distance to the detector (close 
to the front cap to get the highest count rate possible), while the dis- 
tance of the strong 60 Co was varied to create a variable amount of 
background radiation. 



pose LaBr/3:Ce as a well suited candidate material 
to allow y-ray spectroscopy at very high count rates. 
This property has been briefly discussed by comparing 
LaBi^Ce to Nal Q2Q . However, the investigation in [|2|] 
does not include the study of efficiency and resolution 
- the key properties for nuclear physics experiments - 
when increasing the count rate in the LaBr3:Ce detector. 

In this article we investigate in greater detail the high 
count rate capabilities of LaB^Ce. The energy resolu- 
tion as well as the photo peak efficiency are studied with 
respect to the count rate. For further improvement it 
is possible to introduce digital signal-processing (DSP) 
techniques to correct signal distortions and restore ef- 
ficiency losses due to pile-up. It is shown that high- 
resolution y-ray spectroscopy with LaBr3:Ce detectors 
is possible at count rates well in excess of 1 MHz. Sim- 
ulations with synthetic signals have been carried out to 
investigate the limits of this approach, and to support 
the experimental results. 

Section|2]of this article will shortly introduce the ex- 
perimental setup that was used for the measurements, 
while the data analysis, including feature extraction 
from the digitized signals, are detailed in Section [3] 
This section also describes how a timestamp-based pile- 
up correction method, first introduced in [16], can be 
used to further improve the performance of LaBr3:Ce 
detectors at high count rates. In Section [4] the results 
of these investigations will be shown. An outlook for 
future investigations is presented in Section [5] 



2. Detector and Setup 

The experimental setup for the present investigation 
is shown in Fig. Q] The large volume 3"x3" LaBr3:Ce 
detector was irradiated by two different radioactive 
sources. One 137 Cs source of medium activity (45 kBq) 
was placed in contact with the front face of the detector, 
while a 60 Co source with high activity (73 MBq) was 
placed at a variable distance d. The 137 Cs source was 
fixed in place for all measurements and thus provides 



the stable reference, while the strong 60 Co source was 
moved within a range of 5-100 cm in order to vary the 
overall count rate. The count rate using only the 137 Cs 
source was 10 kHz. Bias voltage for the Hamamatsu 
R10233 photomultiplier tube (PMT) was delivered from 
an iseg NHQ-203M power supply through a 2007 volt- 
age divider base from Canberra and operated at bias 
voltages from 300 V to a maximum of 1200 V. The an- 
ode signal of the PMT was amplified by a fast timing 
filter amplifier (TFA) with 1 /is differentiation and 1 ns 
integration time settings in oder to match the dynamic 
range of the sampling ADC. The amplified signal was 
then recorded with a Struck SIS3350 digital sampling 
ADC with a sampling rate of 500 MHz and a resolution 
of 12 bit [17]. The length of each signal trace was set at 
10 5 data points corresponding to 0.2 ms of continuous 
samples. To reach a live time of 2 s, in total 10 4 traces 
have been recorded for each distance d. Data acquisi- 
tion was done using a Struck SIS3104/SIS1 lOOe PCIe- 
VME-interface and a standard PC running the GECKO 
data acquisition software! 18] with online signal pro- 
cessing. Using the PCIe-VME-interface high data trans- 
fer rates of more than 100 MB/s between the VME elec- 
tronics and the PC could be reached to allow fast readout 
of the sampling ADC. 

3. Signal analysis and pile-up correction 

The information conveyed in the digitized signals 
first has to be extracted using DSP methods that are 
described in this section. At first the timestamp is ex- 
tracted from the fast rising edge of the pulse, because 
the steepest part of the slope conveys precise timing in- 
formation. Afterwards regions of the trace without any 
signals are identified based on the timestamps. Then the 
baseline offset is corrected. The energy deposit infor- 
mation can then be extracted from the integral over each 
signal. These parts of the data analysis are described in 
detail in the following. 

3.1. Feature extraction 

Several different steps of signal analysis are per- 
formed on the data in order to extract the basic signal 
information (i.e. energy and timing). The baseline- 
offset voltage is estimated by averaging over 10 data 
points of the signal trace in regions where no y-ray sig- 
nal is present and is then subtracted from the original 
signal. At very high count rates the probability for such 
"empty" regions in a trace decreases, such that the de- 
termination of the baseline becomes more difficult. At 
the cost of precision it is still possible to derive a value 




Figure 2: Top: A differentiating filter (t) transforms the signal 
(dashed, red line) into the timing signal (grey). After filtering with 
an integrating box filter (E), the signal is transformed into a calorime- 
try signal (solid, blue). The delay At between these two signals is 
independent of the amplitude. The trigger time to is determined from 
the timing signal. Bottom: After taking into account the delay At, the 
amplitude of the signal is determined at the time fo. 



for the baseline, if the trace contains at least one of such 
regions. If no suitable regions are present, then differ- 
ent methods have to be employed, which however was 
not necessary in the present case. The dashed red line 
in Fig. [2] shows a typical signal trace with the baseline 
offset voltage already subtracted. To derive a reference 
point in time (trigger time or timestamp) from this sig- 
nal, the data are convoluted with a discrete bipolar Finite 
Impulse Response (FIR) filter kernel with a width that 
equals the signal rise time. It is possible to assume an 
energy independent rise time of about 10 ns, because the 
signal shape does not depend on the energy. This filter 
kernel has the value —1 for the first 10 ns and the value 
+ 1 for another period of 10 ns (see inset "t" in the lower 
panel of Fig. 0. The resulting bipolar timing signal is 
shown in gray. From the maximum of this timing sig- 
nal it is possible to derive an accurate trigger timestamp, 
which is denoted as fo in the illustration. 

In the next step the deposited energy is determined 
by integrating the signal. The original signal is filtered 
with a unipolar FIR filter kernel, i.e. its value is +1 
over its duration, with a width that ideally equals the 
length of the signal (shown with a width of 40 ns in in- 
set "E"), in order to produce an integrated calorimetry 
signal (solid blue). For computational reasons the width 
of the filter can be reduced at the cost of precision. In the 
present analysis a width of 40 ns was used for all of the 
measurements. The amplitude of this calorimetry sig- 



Figure 3: Top: Example of a pile-up situation. Amplitudes (blue) 
have been measured from the calorimetry filtered signal (solid line). 
Bottom: After pile-up correction the measured amplitudes bj are cor- 
rected to yield the true amplitudes a\ (black). The superposition of 
the model pulse shapes (thin lines) results in the filtered signal (thick 
line). 



nal corresponds to the amount of deposited energy. The 
resulting timing and calorimetric signals are aligned in 
time by shifting them with respect to each other, as de- 
noted in Fig.|2]by At. The value of this delay is constant 
for any recorded signal, as long as the signal shape is 
constant, which is the case in the present study. The en- 
ergy deposit is thus given by the value of the calorime- 
try signal at t = to +A t . This approach for extracting the 
signal properties has been chosen for two reasons: First 
this is an exceedingly simple and numerically stable ap- 
proach, and more importantly this approach allows for 
a correction of the measured amplitudes in the case of 
strong pile-up, as described in the following. 

3.2. Pile-Up correction 

At increasing count rates the average time difference 
between two consecutive pulses will be shorter than the 
individual pulse length. Consequently different pulses 
overlap in time, resulting in a correlation of their prop- 
erties, which - if not corrected for - finally leads to a 
reduction of the energy resolution and a loss of photo- 
peak efficiency. The pile-up correction method that was 
used in the analysis has first been presented in II 1611 . The 
detailed description can be found there, and only a short 
outline is given in this article. 

This method relies on the assumption that the pulse 
shape of every signal from the same scintillation detec- 
tor is stable and independent from its amplitude, so that 
the complete signal trace can be expressed as a linear 



combination of single pulses, as illustrated in Fig. [3] 
This assumption is sufficiently well fulfilled. The com- 
plete trace can thus be described by the timestamp of 
each pulse and the corresponding pulse amplitude, as 
well as a single model pulse shape common to all pulses 
that characterizes the dynamics of a single signal. Based 
on this model the correlation of neighboring signals can 
be resolved in order to restore the correct y-ray energy 
information from the measured amplitudes. 

Figure [3] illustrates the principle of the applied pile- 
up correction method. A signal trace (dashed red) with 
five signals from the LaBr3:Ce detector is shown at 
the top together with the computed calorimetry signal 
(blue solid). For illustration purposes the time delay in- 
troduced by the digital filters is removed as described 
above. The single signals partly overlap each other, such 
that the amplitudes b\ taken from the calorimetry signal 
at times f,- signal are overestimated (especially visible 
in the case of b\ and ^5). The solution to correct for 
this overestimation is to treat each of the amplitudes b[ 
as a superposition of the true amplitude a, of the corre- 
sponding signal and additional contributions c, from the 
surrounding signals: 



bi=di + ^Cj. 

=0 



(1) 



j= 



The contribution of the jth neighbor signal can be ex- 
pressed as the true amplitude aj of the neighbor signal 
weighted with a factor my, such that cj = mijcij. The 
value of this factor depends solely on the time differ- 
ence tij = tj — tj between the time of signal i and the 
time of the neighbor j and can be derived with the use 
of the model pulse shape m(t). Taking the value of the 
model pulse shape at the time fy with respect to the trig- 
ger time to of the model pulse yields the desired weight 
rriij = m(to + tij), resulting in 

N N 

bt = £ mj aj = Y< m ( f ° + f y) a .i- ( 2 ) 

j=o ,/=o 

For all signals in the trace this equation can be written 
as a set of linear equations 



b = Ma, 



(3) 



that can be solved by inverting the matrix M containing 
the weights my and thus produce a solution for the true 
amplitudes a. 
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a = M *b 



(4) 
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Figure 4: Top: Generated signal group with a distance of At = 60 ns 
between the two signals. Bottom: Amplitude histograms determined 
for many generated sets of two pulses at At = 60 ns. Distribution of 
first pulse amplitude shown in red (at 0.5), second pulse amplitude in 
blue (at 1.15). 



The model pulse shape is determined using data 
recorded with a low count rate, in order to avoid pile- 
up situations. From this data a large number of single 
signals are extracted. These signals are averaged and 
normalized to amplitude 1.0 to produce the model pulse 
shape. 

The bottom part of Fig. [3] shows the same calorime- 
try signal again, but now together with the computed 
amplitudes a, and the corresponding constituent signals 
(black). These constituent signals are produced by scal- 
ing the model pulse shape with the true amplitude a, as 
determined using Eq.|4]and shifting it in time to match 
the trigger time f,. From this figure it is obvious, that 
the calorimetry signal results from a superposition of 
all constituent signals. The free parameters of this al- 
gorithm are the extents of the matrix M, which control 
how many neighbor signals to take into account for each 
calculation. It is reasonable to assume a small number 
of rows and columns, because the model pulse shape 
only has non-zero elements within about 200 ns around 
its trigger time. This fact allows to calculate the ma- 
trix inversion with only a small amount of computing 
resources, and even enables the implementation of the 
algorithm in hardware for real-time applications (e.g. 
FPGA or DSP as in d). 

3.3. Application limits 

To study the performance of the digital signal analy- 
sis and the limits for the application of the pile-up cor- 



5 1.2 

"ft 1.0 

I 0.8 

St 6 

PS 0.4 



g 1.6 

a i.2 
n 

■a 1-0 

I 0.8 

a 

fi 0.6 

g 0.4 

8 0.2 
Oh 

0.0 



"H 



"Hi! 



"'U, 



II 



■J-f-$-i-i-3 
II 



Hi.: inn.: 

■ I amplitude 

■ 2 amplitude 



■ 1 amplitude 

■ 2'" 1 amplitude 



III 



40 60 80 100 120 140 160 180 200 
At [ns] 



Figure 5: Top: Measured raw amplitudes (&,) with respect to time 
difference At (first signal with squares, second with circles). Hori- 
zontal lines indicate "true" generating amplitudes. Uncertainties in 
each point result from statistical errors. Three distinct ranges in At 
can be identified: (I) Signals are indistinguishable (II) Incorrect am- 
plitudes due to pile-up (III) First amplitude is measured correctly, but 
second one is still overestimated. Bottom: Reconstructed amplitudes 
(fli). Three distinct regions can be identified: (I) Trigger limited mis- 
identification, (II) Transition region with increased uncertainty on the 
result, (III) Incorrectly measured amplitudes are accurately restored. 



rection algorithm, test cases have been constructed with 
two synthetic signals at a well defined distance At. The 
upper part of Fig. [4] shows an example with a distance 
of At = 60 ns and an amplitude ratio of 1:2 of the two 
signals (the second is twice as large as the first). The 
signals are generated from an exponential function with 
a decay constant of 60 ns, which is smoothed with a 
gaussian filter to create a finite rise time of about 20 ns 
(matching the properties of real detector signals) and 
then scaled to its amplitude. White Gaussian noise is 
added to the trace to produce a signal-to-noise ratio 
(SNR) of 10 for the first and 20 for the second signal, re- 
spectively (SNR is here defined as pulse height devided 
by the root mean square of the noise level). This amount 
of noise is chosen to simulate small detector signals. For 
these small signals the limitations of the method with re- 
spect to timing are more visible when compared to large 
signals. The parameter At is varied from to 200 ns in 
steps of 4 ns, and for each step a set of 5000 traces was 
generated. Due to the random noise all generated traces 
in each set differ from each other. The traces were anal- 
ysed with the same approach described in Sec. 13. II to 
determine the trigger times f,- and the raw amplitudes bj 
of the two signals. For each At a histogram (lower part 
of Fig. 3]l was filled with the determined amplitudes, and 



a gaussian distribution was fitted to determine the mean 
value and the FWHM of each of the two peaks. 

The results are shown in the upper part of Fig. [5] as a 
function of the time difference At of the two generated 
signals. The error bars for each point are drawn accord- 
ing to the determined FWHM. Three distinct ranges in 
At can be identified: (I) The time difference is compa- 
rable to or even smaller than the rise time of the signals. 
They appear as one single signal and the trigger algo- 
rithm cannot distinguish between them; (II) At inter- 
mediate time distances, both amplitudes are incorrectly 
measured, because of pile-up. This region extends to 
about three times the rise time for our simulation; (III) 
The first amplitude is measured correctly, because in 
this region the influence of the rising edge of the sec- 
ond signal on the amplitude of the first one becomes 
negligible. The amplitude of the second signal is still 
overestimated (top of Fig. |5j. 

The lower part of Fig. [5] shows the results after the 
pile-up correction method has been applied. For this 
simulation the method only takes into account the next 
neighbors of each signal (l ir order), because only two 
signals participate. The three different regions can be 
identified here as well: Region (I) is the region of unre- 
solvable pile-up, where the pile-up correction cannot be 
applied, because the two signals cannot be separated. 
For the same reason it is not possible to reject these 
signals, based purely on the trigger information. Re- 
gion (II) is the transition region, where the pile-up cor- 
rection method can be applied, and yields the proper 
result. However, for shorter time intervals the statis- 
tical uncertainty increases. Depending on the applica- 
tion, the events in this transition region can be identified 
and rejected. Besides depending on the detector reso- 
lution, these signals may be corrected well enough to 
contribute to the photo peak. The amplitudes in region 
(III) can be properly restored within their intrinsic un- 
certainties. 

The simulation shows that in a range of 60-200 ns 
for At the pile-up correction method yields a good re- 
construction of the signal amplitude, while the uncor- 
rected amplitude of the second signal is determined in- 
correctly. When applying the same method to real data 
an improvement of the extracted properties is expected 
for count rates where a large fraction of the signals have 
a distance in this range. At lower or higher count rates 
the improvement is expected to be much smaller, be- 
cause either an insignificant amount of signals overlap 
or because the average distance between two signals is 
too small to allow for precise determination of the time- 
stamp. 
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Figure 6: (Color online) Energy spectra of emitted y-rays from a 
s0 Co source recorded with the LaBr3:Ce detector. Top: Large dis- 
tance (1.5 m) between source and detector to induce a low count rate 
of 22 kHz. Bottom: The count rate increases to 6 MHz at smaller dis- 
tances (20 cm). Pile-up events are visible above energies of 1.3 MeV 
in the spectrum without correction (dashed red). These are partly re- 
stored in the corrected spectrum (solid blue). 

4. Experimental results 

In this section the experimental results from the mea- 
surements detailed in Sec.|2]are presented. 

Figure[6]shows energy spectra of emitted y-rays from 
the strong 60 Co source recorded with the LaBr3:Ce de- 
tector. For this data only the 60 Co source has been used. 
The top part shows the spectrum taken with a low count 
rate, while the source was placed at a distance of 1.5 m 
to the detector. The lower part of Fig. [6]shows the spec- 
trum when recorded with the source at a distance of 
20 cm. The count rate increases from 20 kHz to 6 MHz 
and pile-up events become apparent in this spectrum, 
which can be seen as high energy tails of the peaks and 
the increase of events with energies much larger than 
1.33 MeV. Without pile-up correction the red (dashed) 
spectrum has less photo peak efficiency (5% peak-to- 
total compared to 16% at 20 kHz) and also the energy 
resolution is worse when compared to the spectrum at 
low count rates (see Fig. Q. After pile-up correction 
taking into account up to three neighboring signals on 
each side (blue) the photo peak efficiency (8% peak-to- 
total) as well as the energy resolution are partly recov- 
ered. At the same time the number of pile-up events 
above 1.33 MeV is reduced to 62%. 

The photo peak efficiency is determined from the 
peak area of the 661 keV peak corresponding to the en- 
ergy of the y-rays emitted by the 137 Cs source. The peak 



Figure 7: Top: Relative resolution at 661 keV ( 137 Cs) without cor- 
rection (red circles) and with correction (blue squares). At the highest 
count rate the resolution without correction could not be determined. 
Bottom: Relative resolution at 1332 keV ( 60 Co). 

region has been fixed for all count rate settings. The ob- 
tained values are then normalized to the peak area at the 
lowest count rate of 28 kHz. 

The dependence of the energy resolution and photo 
peak efficiency on the count rate has been studied quan- 
titatively. The results are presented in Figs. [7] and [8] 
respectively. The energy resolution stays below 4% 
at 661 keV for count rates below 500 kHz, and does 
not increase above 10% for the highest count rate mea- 
sured (14 MHz). The same is true for the resolution at 
1332 keV. From Fig. [7] it is also clear that the pile-up 
correction (blue squares) restores the energy resolution, 
but not by a significant amount. The lengths of the error 
bars correspond to statistical errors only. For the highest 
count rate the uncorrected spectrum could not be anal- 
ysed. 

In the upper part of Fig. [8] the results for the relative 
photo peak efficiency at 661 keV is shown. Again red 
circles indicate the efficiency before pile-up correction, 
and blue squares mark the results after pile-up correc- 
tion. Along with the data points two different curves are 
drawn. These curves correspond to the probability that 
within a certain time interval z a second following sig- 
nal occurs, if a random sequence of signals is assumed. 
These probabilities are calculated using a Poisson dis- 
tribution, according to 



where N denotes the number of signals in the time in- 
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Figure 8: Top: Relative photo peak efficiency at 661 keV normalized 
to the value at the lowest count rate (red circles without correction, 
blue squares with correction). The curves indicate the probability for 
detecting only one signal within a time interval X (dotted: X = 60 ns, 
dashed: X = 200 ns) at the respective count rate. Bottom: Throughput 
(product of relative PP efficiency and count rate). At count rates above 
1 MHz a clear enhancement is visible after correction. The curves are 
calculated from the above probability characteristic. 



terval of length T and n the mean number of pulses per 
time interval (count rate). Thus, P(0,n) is the prob- 
ability for a second pulse occuring within the time T 
for a given count rate n. For the illustration two dif- 
ferent values for T (blue: T = 60 ns, red: T = 200 ns) 
as estimated from the simulation in Sec. 13.31 are used. 
The measured data follow the shape of the two curves 
in both cases, which indicates that the simulation cor- 
rectly predicts the region where the application of the 
pile-up correction method will improve the results. In 
the range between 500 kHz and 10 MHz the advantage 
from using the pile-up correction method is most no- 
ticeable. For lower count rates the effect of pile-up can 
be neglected, while for higher count rates the leading 
edges of the pulses overlap and thus prohibit accurate 
triggering. The lower part of Fig.[8]shows the resulting 
estimated throughput (given by the product of relative 
efficiency and the count rate) together with the ideal 
case where the throughput equals the incoming count 
rate (black line) and the two curves from the theoretical 
prediction. This figure demonstrates, that at a low count 
rate the effect of pile-up correction is negligible, while 
at count rates above 1 MHz the improvement is signif- 
icant, increasing the ideal count rate and the maximum 
throughput by a factor of 3. 



5. Conclusions and Outlook 

The analysis has shown that LaB^Ce detectors are 
suitable for y-ray spectroscopy even at very high count 
rates, exceeding 1 MHz. Using the presented pile-up 
correction method, the loss of photo peak efficiency 
from piled-up pulses can be alleviated. A reasonable 
range for the application of this method is between 
500 kHz and 10 MHz, corresponding to an average dis- 
tance of two signals in the range from 60 to 200 ns. At 
count rates above 10 MHz the analysis of the digitized 
signals becomes more demanding, because additional 
effects decrease the accuracy of feature extraction. One 
of these effects is the saturation of the PMT with in- 
creasing count rate, which can be avoided either with 
a lower bias voltage or by using an active voltage di- 
vider with transistor stages. Other effects that arise at 
these high count rates like unresolvable pile-up and the 
fact that the baseline offset can not be determined in a 
simple manner need additional investigation. 
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